Cancer Cell
○ Elsevier BV
Preprints posted in the last 7 days, ranked by how well they match Cancer Cell's content profile, based on 42 papers previously published here. The average preprint has a 0.05% match score for this journal, so anything above that is already an above-average fit.
Lam, J. M.; Walker-Samuel, S.; Pennycuick, A.
Show abstract
Somatic copy-number amplification is pervasive in cancer, and the genes it carries are candidate drug targets - but only those whose amplification is transmitted to accessible surface protein can be reached by an antibody-drug conjugate (ADC). We build an integrated map of copy-number-to-protein transmission across six tumour types and ask, for every amplified gene, whether its dosage reaches the surface. Copy number transmits to mRNA (median per-gene r = 0.21) but is attenuated at the protein level in 85% of genes, and the mRNA ranking is largely preserved to protein (rho = 0.70); the ranking is set principally at the chromatin/transcription step - among directly measured regulatory inputs, promoter DNA methylation and tumour chromatin accessibility each explain about an order of magnitude more of the transmission variance than gene structure, and do so complementarily. Critically, transmissibility is a stable, gene-intrinsic property: it is predictable from gene properties alone, with no proteomic input, at a leave-gene-out rank correlation of 0.52 (R2 = 0.29); it is not positional (holding out whole chromosome arms changes accuracy by 0.001); and it transfers across lineages (Kendall W = 0.97 across leave-one-lineage-out refits). This licenses a predictor that nominates surface targets in cancer types that lack a tissue-referenced proteome, combining direct protein measurement where it is available with prediction where it is not. Requiring co-elevation on a recurrent amplicon with measured transmissibility and an accessible extracellular ectodomain nominates 22 surface antigens on 18 distinct recurrent amplicons across four cancer types (renal, endometrial and both lung subtypes) - for example ITGB8+TSPAN13+TTYH3 on lung 7p, NCSTN+HSD17B7+MPZL1 on 1q (recurrent in several types), the transferrin receptor TFRC on squamous 3q, and FZD1 on clear-cell renal 7q; 21 of the 22 are non-driver passengers and 10 are confirmed on the experimental Cell Surface Protein Atlas. In single malignant cells, against a null that controls for per-cell sequencing depth, the co-detected constructs sit at a modest 1.05-1.45x above independence (p < 0.001, donor-block bootstrap intervals clear of 1.0), and at binding-relevant thresholds the normal-tissue co-expression collapses - so an avidity AND-gate that binds stably only where the antigens co-occur would spare normal cells that carry only one. Observed transmissibility itself transfers strongly between the two lung subtypes ({rho} = 0.88) and remains positive across distant lineages, consistent with the shared cell-of-origin regulation the map implies. Single-cell co-detection is demonstrated wherever a malignant single-cell atlas exists (both lung subtypes and glioblastoma - the latter entirely from prediction, using no GBM surface-abundance measurement); the remaining cohorts are nominated on the same genetic and topological evidence. The result is a pan-cancer, confidence-tiered catalogue of multi-antigen ADC co-target sets with a concrete plan to test them.
Uppalapati, S. C.; Butler, D. W.; Bouobda, G.; Liptrap, E. J.; Schmalz, P. G.; Holland, M. T.; Riley, K.; Filippova, N.; Nabors, L. B.; Markert, J. M.
Show abstract
Background: Glioblastoma remains resistant to most immune-based therapies. Surgery may create a perioperative window in which systemic immune activation and tumor antigen release intersect. We evaluated whether COVID-19 vaccination shortly before first glioblastoma surgery was associated with survival. Methods: We performed a retrospective single-center cohort study of adults with newly diagnosed glioblastoma undergoing initial biopsy or resection from 2021 to 2025. The primary exposure was documented COVID-19 vaccination within 100 days before first tumor surgery. Overall survival was analyzed from surgery using Kaplan-Meier and Cox models, with 1:1 propensity matching and sensitivity analyses addressing treatment completion, calendar time, surgical selection, steroid exposure, immune-cell variables, COVID severity, and negative-control vaccination. Results: The cohort included 187 patients: 64 perioperatively vaccinated and 123 non-perioperative comparators. Among vaccinated patients, 59/64 (92.2%) received mRNA vaccines; median vaccination-to-surgery interval was 81 days (IQR 71-90). Median overall survival was 743 days in vaccinated patients versus 318 days in comparators (unmatched HR 0.48, 95% CI 0.30-0.76; p=0.002). After 1:1 matching, median survival was 743 versus 349 days (HR 0.52, 95% CI 0.34-0.80). Sensitivity analyses accounting for adjuvant therapy, surgery year, extent of resection, steroid exposure, immune-cell measures, and COVID hospitalization were directionally consistent. Influenza vaccination was not associated with survival. Conclusions: COVID-19 vaccination within 100 days before first glioblastoma surgery was associated with longer overall survival. These findings identify perioperative vaccination timing as a potentially relevant and modifiable variable in glioblastoma outcomes.
Manjarrez, S.; Diaz, F. C.; Carranza, F. G.; Waldrup, B.; Ninova, M.; Velazquez-Villarreal, E.
Show abstract
Background: Early-onset colorectal cancer (EOCRC) is increasing globally, particularly among Hispanic/Latino (H/L) populations, yet the contribution of tumor-colonizing microbiota to age-associated colorectal cancer (CRC) biology remains poorly understood. Most microbiome studies have focused on fecal communities or non-Hispanic populations, leaving the intratumoral microbial landscape of H/L patients largely unexplored. Methods: We performed an exploratory characterization of tumor-colonizing microbiota using whole-exome sequencing (WES) data from four primary colorectal tumors obtained from H/L patients treated at City of Hope, including two EOCRC (<50 years) and two late-onset colorectal cancer (LOCRC; [≥]50 years) cases. Following removal of host-derived sequences, microbial taxonomic profiling was conducted at the family, genus, and species levels, and microbial metabolic pathways were inferred. Clinical and pathological data were integrated to evaluate age-associated differences in microbial composition and predicted function. Results: Family-, genus-, and species-level analyses consistently demonstrated greater microbial diversity in LOCRC than EOCRC. LOCRC contained more than twice the number of unique bacterial families, nearly three times as many unique genera, and more than twice as many unique bacterial species. A conserved core microbiota, including Fusobacteriaceae, Prevotellaceae, Fusobacterium, and Prevotella, was identified across both age groups, whereas LOCRC was enriched in CRC-associated taxa including Fusobacterium nucleatum, Bacteroides fragilis, Parvimonas micra, Porphyromonas asaccharolytica, and Dialister pneumosintes. Species-level analyses revealed only a single shared bacterial species between EOCRC and LOCRC, indicating progressive microbial divergence with increasing taxonomic resolution. In contrast, functional profiling identified 11 predicted microbial metabolic pathways, of which nine were shared between age groups, two were unique to EOCRC, and none were exclusive to LOCRC. Core metabolic pathways involved in energy metabolism, amino acid biosynthesis, phospholipid metabolism, and central carbon metabolism exhibited comparable abundance across both groups, demonstrating substantial functional conservation despite pronounced taxonomic differences. Conclusions: Tumor-colonizing microbiota differ markedly between EOCRC and LOCRC in H/L patients, with late-onset tumors exhibiting substantially greater microbial richness and taxonomic complexity. Despite these compositional differences, microbial metabolic functions remain largely conserved, supporting the concept of functional redundancy within the colorectal tumor microenvironment (TME). Although exploratory, this proof-of-concept study provides one of the first characterizations of intratumoral microbiota in H/L EOCRC and establishes a foundation for larger multi-omics investigations aimed at identifying microbiome-based biomarkers and therapeutic targets for precision oncology.
Sanfeliu, E.; Segui, E.; Martinez-Romero, A.; Albarran-Fernandez, V.; Pascual, T.; Marin, M.; Martinez-Saez, O.; Gomez-Bravo, R.; Garcia-Fructuoso, I.; Rodriguez-Hernandez, A.; Walbaum, B.; Galvan, P.; Angelats, L.; Rubio-Perez, C.; Saura, C.; Oliveira, M.; Ciruelos, E.; Manso, L.; Pernas, S.; Vidal, M.; Waks, A. G.; Tolaney, S. M.; Pare, L.; Parker, J. S.; Villagrasa, P.; Ferrero-Cafiero, J. M.; Perou, C. M.; Campo, E.; Tabernero, J.; Braso-Maristany, F.; Prat, A.
Show abstract
Tumor-infiltrating lymphocytes (TILs) are widely used to assess antitumor immunity in breast cancer but may not reflect the functional competence of adaptive immune responses. We show that immune organization, reflected by tertiary lymphoid structures (TLS) and coordinated humoral and cellular immune programs, represents a distinct dimension of tumor immunity beyond lymphocyte abundance. By integrating histologic, transcriptomic, spatial, and immune receptor profiling analyses across multiple breast cancer cohorts, we show that immune organization is associated with greater immune repertoire diversity, evidence of therapy-induced clonal selection, and improved clinical outcomes, independent of immune infiltration. Transcriptomic measures of immune organization retained independent prognostic value across external cohorts, whereas measures of immune infiltration did not. Furthermore, treatment-induced increases in immune organization, but not immune infiltration, were associated with therapeutic response. These findings identify immune organization as a dynamic and clinically measurable state of adaptive antitumor immunity with implications for prognosis, treatment monitoring, and therapeutic development in breast cancer.
Rentroia-Pacheco, B.; Sharma, H.; Pozza, L.; Traets, J. J. H.; Tandukar, B.; Steijlen, O. F. M.; Ruiter, R.; Cruz-Pacheco, N.; Huigh, D.; Van Hoeck, A.; Chen, Y.-T.; Infante, B.; Baskurt, D.; Arunachalam, V.; Eggermont, C. J.; Bas-Cristobal Menendez, A.; Nijsten, T.; van de Werken, H. J. G.; Mooyaart, A. L.; Bellomo, D.; Wakkee, M.; Shain, A. H.; Hollestein, L. M.
Show abstract
Cutaneous squamous cell carcinoma (cSCC) is the second most common form of cancer worldwide. While most cSCCs are not life-threatening, 2-5% of patients develop metastases. To better understand what causes some cSCCs to progress to metastatic disease, we assembled a nationwide cohort of 19,120 patients with clinico-pathologically annotated tumors linked to metastatic outcome. RNA-sequencing was performed on 378 tumors, and whole-exome sequencing on 147, with balanced numbers of tumors that progressed to metastatic disease (cases) and did not (controls). UV radiation was the dominant mutational signature with additional contributions from aging, APOBEC activity, and, in immunosuppressed patients, azathioprine exposure. We identified 38 genes under selection across a core set of signaling pathways. Gene expression clusters were primarily associated with the differentiation state of tumor cells and secondarily with the composition of the tumor microenvironment. Several mutational and transcriptional programs were associated with metastasis, including a dedifferentiated gene expression signature, activating mutations in the RAS signaling pathway, loss-of-function alterations in the SWI/SNF chromatin remodeling complex, and specific arm-level copy number alterations. A 23-gene expression signature was built to predict metastasis from primary cSCC tissue. The signature was validated in two independent cohorts (N=102 and 52), where it predicted metastasis independently of staging systems. Together, these findings provide the most detailed molecular portrait of cSCC to date and establish an assay for risk stratification suitable for clinical implementation.
Su, Z.; Li, T.
Show abstract
The therapeutic landscape for hepatocellular carcinoma (HCC) is evolving rapidly, necessitating scalable approaches to synthesize the expanding scientific literature. We characterized thematic shifts in HCC treatment and prognosis research by conducting a retrospective bibliometric analysis of influential publications from 2023 and 2024. Using the OpenAlex database, we identified the 50 most highly cited papers from each year based on eighteen-month post-publication citation counts. Large language models were deployed to extract, normalize, and classify concepts from unstructured text into canonical topics and parent themes, enabling quantitative year-over-year frequency comparisons. Analysis of these 100 papers revealed a distinct maturation in research focus. Although broad categories like general immunotherapy remained prevalent, their relative frequency declined in favor of specific dual immune checkpoint regimens, notably CTLA-4 inhibition and the durvalumab plus tremelimumab combination. Concurrently, parent themes related to radiomics, imaging, and health systems exhibited significant growth in the 2024 cohort. These findings demonstrate a thematic transition in high-impact HCC research from foundational immuno-oncology toward optimized combination therapies and precision diagnostics. Furthermore, this study highlights the utility of artificial intelligence-driven bibliometrics for objectively tracking dynamic conceptual shifts in oncology. A web interface for exploring the data is available at https://pri.pepkio.com/.
Kumar Reddy, K.; Hahn, W.; Winter, S.; Roellig, C.; Mueller-Tidow, C.; Serve, H.; Baldus, C. D.; Fransecky, L.; Schliemann, C.; Burchert, A.; Schaefer-Eckart, K.; Kaufmann, M.; Schetelig, J.; Bornhaeuser, M.; Middeke, J. M.; Eckardt, J.-N.
Show abstract
Rising costs, slow accrual and molecular substratification of cancers necessitate novel clinical trial designs. We demonstrate that artificial intelligence-generated synthetic patients can replace real controls to reproduce results of the SORAML trial. Using external multimodal data from 1,377 acute myeloid leukemia (AML) patients from previous trials and a real-world registry, we fine-tuned a tabular foundation model to generate synthetic patients, reproducing clinical and genetic features and outcome associations. Synthetic patients were then matched to the original SORAML intervention group using Cox risk scores, replacing the original control and reproducing the original trial result with near-identical median event-free survival (EFS) and treatment effect (original hazard ratio [HR] 0.64, 95%-confidence interval [CI] 0.47-0.87, p=0.004; with synthetic control HR 0.66, 95%-CI 0.48-0.90, p=0.009). Our findings demonstrate that AI-generated synthetic patients can serve as statistically rigorous controls supporting novel trial designs.
Collins, M. P.; Lahr, D. L.; Topal, S.; Khalil, A.; Hickman, D.; Spidale, N.; Pandit, N.; Reilly, S.; Lyons, K.; Horrigan, K.; Zhao, T.; Batonga, J.; Bosinger, M.; D'Aco, K.; Ball, B.; Kishtagari, A.; DiNardo, C. D.; Stein, E. M.; Quintas-Cardama, A.; Smolen, G. A.
Show abstract
Impaired cellular differentiation is a defining characteristic of myeloid malignancies and remains a major therapeutic challenge. The BRG1/Brahma-associated factor (BAF) chromatin remodeling complex, through the ATPases SMARCA4 and SMARCA2, maintains the stemness of leukemic blasts and thus represents a promising target for novel differentiation-based therapies. In a phase 1 study in advanced myeloid malignancies, the first-in-class dual SMARCA4/2 inhibitor FHD-286 combined with decitabine (DAC) was tolerated and produced an objective response rate of 12.8% (6/47) compared with no responses with FHD-286 monotherapy. To understand the basis of this activity, we integrated high-dimensional flow cytometry and single-cell genomic analyses of longitudinal bone marrow samples from responders and nonresponders. While FHD-286 monotherapy was predominantly associated with myeloid differentiation, responders to FHD-286+DAC combination therapy exhibited a range of myeloid and erythroid differentiation trajectories. FHD-286 potentiated the transcriptional impact of DAC, driving tumor clones to fully differentiate out of the immunophenotypically and transcriptionally defined blast compartment. Responders had a baseline transcriptional profile similar to that of CEBPA-mutant acute myeloid leukemia and showed further downregulation of CEBPA upon treatment. These findings reinforce tumor cell differentiation as a mechanism of response to pharmacologic SMARCA4/2 inhibition and support further evaluation of FHD-286+DAC in molecularly defined patient subsets.
Alford-Holloway, M. N.; Reed, S. C.; Pershad, Y.; Van Amburg, J. C.; Potts, C.; Mohan, S. R.; Luo, L. Y.; Ferrell, P. B.; Savona, M. R.; Park, B. H.; Johnson, D. B.; Bick, A. G.; Kishtagari, A.
Show abstract
Background The clinical significance of clonal hematopoiesis of indeterminate potential (CHIP) in melanoma remains incompletely defined, particularly with respect to CHIP genotype, clone size, and somatic mutations (e.g BRAF mutations). We integrated human cohort data and a syngeneic melanoma mouse model to evaluate whether CHIP is associated with melanoma risk, tumor growth, and differential clinical outcomes. Methods We analyzed CHIP prevalence and survival in a large treatment-unselected melanoma cohort (n=2,480), evaluated tumor growth in a syngeneic BRAF-mutant (BRAFmut) melanoma murine model of TET2-CHIP and DNMT3A-CHIP, and assessed survival outcomes in an immune checkpoint inhibitor (ICI)-treated advanced melanoma cohort (n=361). Associations with progression-free survival (PFS) and overall survival (OS) were evaluated using Kaplan-Meier analyses and multivariable Cox proportional hazards models. Results CHIP was enriched among patients with treatment-unselected melanoma compared with age/sex-matched healthy controls, and larger CHIP clone size showed an age-adjusted association with inferior OS. In a syngeneic BRAFmut melanoma murine model, TET2-CHIP, but not DNMT3A-CHIP, was associated with significantly increased primary melanoma tumor growth. Among patients with ICI-treated advanced melanoma, CHIP was associated with worse OS compared with patients without CHIP. TET2-CHIP had the strongest adverse association with survival, whereas DNMT3A-CHIP was not significantly associated with PFS or OS. Conclusions CHIP is enriched in melanoma and exploratory analyses demonstrate genotype-specific differences in melanoma tumor growth and clinical outcomes. These findings support further investigation of genotype-specific CHIP profiling as a potential biomarker for melanoma risk stratification and immunotherapy outcomes.
Slotman, E.; van Disseldorp, L. M.; de Jong, G.; Fransen, H. P.; Reyners, A. K. L.; Tol, J.; Jager, A.; Westgeest, H. M.; Sonke, G. S.; van Laarhoven, H. W. M.; van Zuylen, L.; van den Heuvel, M. M.; Koopman, M.; Smit, E.; Raijmakers, N. J. H.; Siesling, S.
Show abstract
Introduction: This study aimed to provide population level survival trends during the era in which new systemic therapies transformed treatment guidelines for metastatic cancer, as well as insights on the real world use of these treatments and associated survival. Methods: Adults diagnosed with synchronous metastatic solid cancer in 2008 until 2022 (22 cancer types) were identified from the Netherlands Cancer Registry. Median overall survival (OS) was assessed by five year diagnostic period. For 2018 until 2022, systemic therapy use in any treatment line was analyzed and survival percentiles within treatment and cancer types were estimated with Kaplan Meier survival analyses. Results: Median OS in the overall cohort (n=280,419 patients) improved from 6 to 8 months between the period 2008 until 2012 and 2018 until 2022. Among patients diagnosed in 2018 until 2022, 15% received immunotherapy, 15% targeted therapy, 29% chemotherapy and/or traditional hormone therapy only, and 39% no systemic therapy. In some cancer types, a relatively large proportion of treated patients had longterm survival (e.g., immunotherapy in melanoma: p50 = 67 months). Other cancer types had a smaller subset of treated patients (p10 and p25) with substantially better outcomes than the median (e.g., targeted therapy in NSCLC: p50 = 22 months, p10 = 96 months). Conclusion: Population level survival for patients with synchronous metastatic solid cancer has modestly improved over time. The marked survival heterogeneity within cancer and treatment types highlights both the potential and uncertainty associated with (novel) treatments. Improved prediction of treatment effects and clear communication regarding survival expectation remain critical. Presenting multiple survival scenarios over median survival alone can support decision making.
Bielcikova, Z.; Tichopad, A.; Rybar, M.; Petrakova, K.; Rozanek, M.; Mothejlova, K.; Dusek, L.; Donin, G.
Show abstract
Population-based mammography screening improves breast cancer outcomes, but its impact on real-world treatment pathways and quality indicators (QIs) remains incompletely described. We conducted a retrospective nationwide cohort study using linked data from the Czech National Cancer Registry and the National Registry of Reimbursed Health Services. Women aged [≥]18 years with a first breast cancer diagnosis between 2017 and 2024 were classified as screen-detected (SCR) or diagnostically-detected (DIG) according to the imaging modality preceding histological verification. Outcomes included stage distribution, untreated cases, first-line treatment, main treatment modality, time to treatment, multidisciplinary team discussion (MDT), centralization to Comprehensive Cancer Centres (COCs), and survival patterns. The verified cohort included 47,648 women: 26,817 SCR cases (56.3 %) and 20,831 DIG cases (43.7 %). In this nationwide analysis, SCR breast cancer was associated with earlier stage at diagnosis and better survival patterns, but also with longer time to treatment and longer time to MDT discussion than DIG-detected disease. Although treatment rates were high and centralization improved over time, substantial regional variation persisted in care pathways, MDT use, and access to COCs. These findings support continued strengthening of screening participation, monitoring of care intervals, and quality assurance of MDT reporting and regional oncology care delivery.
Madrigal, A.; Kim, M.; Mehrjoo, Z.; Nishimura, T.; Saatci, O.; Osakwe, A.; Zavacky, E.; Moslemi, E.; Glennon, K. I.; Dankner, M.; Maritan, S. M.; Kuasne, H.; Pilon, V.; Monast, A.; Soytas, M.; Arseneault, M.; Oikonomopoulos, S.; Harutyunyan, A.; Lu, T.; Rayes, R.; Soto, L. M.; Hernandez-Corchado, A.; Spicer, J. D.; Petrecca, K.; Siegel, P.; Park, M.; Ragoussis, J.; Sahin, O.; Brimo, F.; Tanguay, S.; Riazalhosseini, Y.; Najafabadi, H. S.
Show abstract
While extensive cellular heterogeneity in renal cell carcinomas (RCC) is linked to diverse clinical outcomes, our understanding of this diversity is limited to those driven by clonal patterns or activity of canonical pathways. Here, we present a compendium of over 85,000 single-cell gene expression profiles from primary and metastatic tumors as well as patient-derived models across four RCC subtypes, including the rare clear cell papillary renal cell tumors, which we show are often misclassified and for which we identify CASP14 as a highly sensitive and specific biomarker. We dissect malignant cell variation within and across tumors using a generative modeling framework that accounts for clonal and copy number-driven expression shifts, defining 59 gene expression programs that deconstruct canonical pathways into functional submodules with divergent activity patterns, distinct regulators, and differential association with clinical outcomes. Despite the canonical view that VHL-deficient clear cell RCC exists in a constitutive pseudohypoxic state, we show strong intra-tumor variability of a hypoxia inducible factor 2 (HIF2)-driven program linked to poor outcome. We also identify early, spatially organized activation of a complete epithelial-to-mesenchymal transition (EMT) program, loss of epithelial identity, and upregulation of protein translation programs as key characteristics of metastatic progression. Finally, a metastatic signature capturing cellular de-differentiation and translational activity identifies primary tumors associated with adverse clinical outcomes. Together, this resource establishes a framework for dissecting malignant cell heterogeneity, refines RCC subtype classification, and defines transcriptional programs underlying metastasis progression.
Kamelian, K.; Pascall, D. J.; Cheng, M. T. K.; Meng, B.; Altaf, M.; Morse, R. M.; Aggio, J. B.; Egan, D. J. S.; Chen-Xu, M.; Trivioli, G.; Sutton, B.; Richter, A.; Gonzalez-Vazquez, L. D.; Cormie, C.; Kemp, S.; Yeadon, R.; Hyatt, B.; Wong, A.; Thesin Pelamkulangara, N.; Fraser, E.; McCarthy, B.; Novaes, F.; Stott, S.; Galvin, A.; Bellis, K. L.; De Angelis, D.; Harrison, E. M.; Martin, D.; Smith, R. M.; Gupta, R. K.
Show abstract
Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.
Weerasinghe, C.; Osowicki, J.; Simpson, J. A.; Crocker-Buque, T.; McCarthy, J.; Williams, E.; Price, D. J.
Show abstract
Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.
Prosty, C.; Butler-Laporte, G.; Brophy, J.; Frenette, C.; Loo, V.; Coburn, B.; Hota, S.; Longtin, Y.; Kong, L.; Muller, M.; Steiner, T.; Valiquette, L.; Daneman, N.; Daley, P.; Nott, C.; MacFadden, D. R.; Kandel, C.; Chen, Y.; Perez- Patrigeon, S.; Lee, T. C.; McDonald, E.
Show abstract
Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.
Brochu, H. N.; Shi, Q.; Song, K.; Zhang, Q.; Munroe, J.; Harris, N. J.; Britt, N.; Zeng, Q.; Kapuria, K.; Chappell, J.; Norvell, B. M.; Peavy, L.; Williams, J. D.; Harris, A. B.; Chaitram, J.; Hutson, C. L.; Deng, J.; McGrath, D.; Boles, D.; Dale, S. E.; Gigante, C. M.; Iyer, L. K.
Show abstract
Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.
Nimalrathna, S. U.; Harischandra, H.; Kimber, M.; Chandrasena, N.; De Silva, N.; Mallawarachchi, H.; De Silva, B. G. D. N. K.
Show abstract
The World Health Organization (WHO) validated Sri Lanka had eliminated lymphatic filariasis as a public health problem in 2016, the second country in Southeast Asia to attain this status. However, post-validation surveillance has identified sporadic cases of brugian filariasis. The reemergence of Brugia malayi infections in Sri Lanka warrants urgent investigations. Recent studies have shown that the parasite responsible for the reemergence is a novel zoonotic Brugia sp. maintained among dogs that is closely related but distinct to the human-infecting B. malayi species. The current study employed morphological and morphometric assessments, revealing that this novel zoonotic Brugia sp. is within the B. malayi morphological range. Molecular characterization of three genomic regions, the nuclear genomic region SLXI, the non-coding region HhaI, and the mitochondrial genomic region COXI confirmed it as a genetic variant more closely related to B. malayi than to B. pahangi. Phylogenetic analysis further indicated it as a distinct genomic variant, closely related to a B. malayi-like parasite reported from India. Notably, that same parasite was identified in infected humans, animals, and potential vector mosquitoes. This, together with the detection of both human and animal blood within the same brugian infective mosquitoes, and delineating the canine origin of the parasites in human infections, provides compelling evidence supporting zoonotic transmission of this parasite. To our knowledge, this is the first report demonstrating the presence of the same brugian parasite in humans, domestic animals, and potentially infective mosquitoes in Sri Lanka, supported by multi-genomic evidence. The recent identification of multiple potential mosquito vector species suggests that this parasite may have undergone adaptive changes, facilitating its ability to overcome the species barrier. These findings substantiate the long-held hypothesis of zoonotic transmission of the reemerged brugian parasite, highlighting significant implications for ongoing surveillance and control strategies.
Gu, S.; Petrovitch, D.; Hall, O. T.; Lambert, J. W.; Kember, R. L.; Nahid, N. A.; Ma, Q.; Sprague, J. E.; McDonough, C. W.; Johnson, J. A.
Show abstract
Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.
Liu, J. B.; Chen, Y.-J.; Edelen, M. O.; Pusic, A. L.; Martin, N. E.; Zeng, C.
Show abstract
Purpose: Nonresponse to routinely collected patient-reported outcome measures (PROMs) threatens the representativeness of aggregated data. We characterized patient-, provider-, and clinic-level factors associated with PROMIS Global-10 nonresponse in routine radiation oncology care. Methods: In this retrospective cohort study, all adults seen at five Mass General Brigham radiation oncology clinics over one year were included. The primary outcome was patient-level nonresponse, defined as never completing the portal-administered Global-10 versus completing it at least once. Using iterative mixed-effects logistic regression, we modeled patient-, provider-, and clinic-level factors. Results: Among 12,214 patients, 71 providers, and five clinics, patient- and appointment-level response rates were 35.4% and 10.9%, with patient-level response ranging nearly fivefold across clinics (12.8% to 66.2%). In Model 1, male sex, lower education, not working, and recent surgery had higher odds of nonresponse, and longer time since diagnosis lower odds. After provider- and clinic-level factors were added, patient sex, education, and employment became nonsignificant, whereas recent surgery (adjusted odds ratio [aOR] 1.97) and longer time since diagnosis (aOR 0.46 for >12 months) persisted. A provider's historical collection rate was protective but attenuated at the clinic level. There, a later program launch (aOR 0.29) and higher historical collection rate (aOR 0.79) correlated with lower nonresponse, whereas academic versus community setting did not. Conclusions: Nonresponse to routinely collected PROMs is a multilevel phenomenon driven substantially by clinic-level implementation factors, not patient characteristics alone. Because response rate is only a proxy for representativeness, PROMs programs and PRO-based performance measures should prioritize representative collection over volume.
Djaafara, B. A.; Elyazar, I. R.; Yosephine, P.; Surya, A.; Silalahi, F. S.; Handito, A.; Thohir, B.; Aryani, D.; Gunawan, D.; Nisa, A. K.; Prianto, E.; Samad, I.; Cook, A. R.; Huang, A. T.; Clapham, H. E.; Bhatt, S.; Mishra, S.
Show abstract
Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.